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of fractional derivative type 

Teodor M. Atanackovic*, Stevan Pilipovic^ and Dusan Zorica* 

July 3, 2012 

Abstract 

We study forced oscillations of a rod with a body attached to its free end so that the motion 
of a system is described by two sets of equations, one of integer and the other of the fractional 
order. To the constitutive equation we associate a single function of complex variable that 
plays a key role in finding the solution of the system and in determining its properties. This 
function could be defined for a linear viscoelastic bodies of integer/fractional derivative type. 

<- ' Keywords: fractional derivative, distributed-order fractional derivative, fractional vis- 

Q j coelastic material, forced oscillations of a rod, forced oscillations of a body 

1 Introduction 

In this paper we continue our recent work on the dynamics of viscoelastic rods described through 
the fractional derivative type equations, presented in [H [SJ El [JJ [5] ■ A problem that we shall be 
dealing with in the present paper is forced oscillations of a body attached to a viscoelastic rod 
with comparable masses. A rod-body system is shown in Figure [I] 
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jj^ ■ Figure 1: System rod-body. 

Our new approach in the investigation of the dynamics of linear viscoelastic rods of fractional 
type, proposed in [5], is based on the properties of a specially denned function of complex 
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variable M (see (fT5| ). Function M is associated with the Laplace transform of the constitutive 
equation for the material of a rod. It is defined similarly as complex moduli (a quantity obtained 
after application of Fourier transform to constitutive equation). By considering two cases of 
constitutive equations, we show in this paper (that is a continuation of [8]) that the proposed 
approach can be adventitiously used in dynamics of viscoelastic rods. 

We find an explicit form of the solution as a convolution of a forcing function and a solution 
kernel. Moreover, we present numerical examples, corresponding to two common cases of con- 
stitutive equations. We analyze forced oscillations of a system consisting of a viscoelastic rod 
of fractional-order type and a body attached to its end. Thus, the cases of the dynamics of an 
elastic rod and of a light rod (mass of a rod is negligible) are special cases of our analysis. We 
refer to [TS] for the analysis of oscillations of an elastic rod with the mass attached to its end. 

In [3], we analyzed forced oscillations of a body attached to a viscoelastic rod, described by a 
fractional distributed-order model. We assumed that the mass of a rod is negligible compared to 
the mass of a body. Similarly as in [3] , we analyzed in [SJ [7] the wave propagation in a viscoelastic 
solid-like rod of finite length with one of its ends fixed to a rigid wall. We considered two cases 
(i.e. two types of boundary conditions): the case when there is a prescribed displacement and 
the case when there is prescribed stress on rod's free end. Similar problem of wave propagation 
was analyzed in [5] for a rod made of viscoelastic fluid-like material. The present paper is closely 
related to j] where we analyze a more general form of a constitutive equation. Actually, in [5] we 
gave a theoretical background which we use here and analyze two models which will be described 
below. 

Let m be the mass of a body attached to a rod. The length of the rod in undeformed state 
is L and its axis, at the initial time moment as well as during the motion, coincides with the 
x axis, see Figure [TJ Let x denote a position of a material point of the rod at the initial time 
to = 0. The position of this point at the time t > is x + u (x, t) . The equations of motion of 
the rod-body system are 

d d 2 d 

—a(x,t) = p- K p?u(x,t), e (x, t) = — u (x, t) , x £ [0, L] , t > 0, (1) 
ox ot z ox 

1 ,1 

CT ( 7 ) B]a(x,t)dj = E / e ( 7 ) D7e(a: > i)d7, i6[fl,L], t > 0, (2) 
o Jo 
d 

u(x,0)=0, —u(x,0)=0, a(x,0)=0, e(x,0)=0, i6[0,I], (3) 

u(0,t)=0, -Aa(L,t)+F(t) =m—u(L,t), t > 0. (4) 

ot 

We use symbols a, u and e, in the equation of motion ([TJ) x and in the strain ([TJh, to denote 
stress, displacement and strain, respectively, depending on the initial position x at time t, while 
p denotes the density of a material. Constitutive equation ([2]) corresponds to the distributed- 
order fractional derivative model of a viscoelastic body, E is a generalized Young modulus (a 
positive constant having the dimension of a stress), 4> a and <j) E are given constitutive functions or 
distributions, oDjf is the left Riemann-Liouville fractional derivative operator of order 7 € (0, 1) 
(see EI]) 

oD7y(*):=^(r^)*y(*)), *>o, 

where Y is the Euler gamma function and * is the convolution. Recall, if /, g G Lj oc (R) , 
supp /, g C [0, 00) , then (/ * g) (t) := f* f (r) g (t - t) dr, tel. We refer to QS QS E] for 
the basic definitions and assertions of fractional calculus. Initial conditions in ([3]) specify that 
the rod-body system is unstressed and in the state of rest at the initial time instant. Boundary 
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condition ((][]) i means that one end of the rod is fixed. The other boundary condition ^2 is the 
equation of translatory motion along the x axis of the body attached to the free end of the rod. 
In (H]) , A stands for the cross-section area of the rod and F stands for the known external force 
acting on the body. 

Regarding the constitutive equation we consider the following two cases. 

I Fractional Zener model of a viscoelastic body 

(1 + a D«) a (x, t) = E (1 + b Df) e (x, t) . (5) 

It is obtained from ^ by choosing 

<l>Al)~Hl) + aS(j-a), e ( 7 ):=*( 7 ) + 6<y(7-a) ) a e (0, 1) , < a < 6, (6) 

where S denotes the Dirac delta distribution. 

II Distributed-order model of a viscoelastic body 

a" / D]a(x,t)d-f = E [ b 1 ^Dje (a, t) d 7 . (7) 
Jo 

It is obtained from @ by choosing 

ff (7) := a\ <j> £ (7) := b\ 7 G (0, 1) , < a < b. (8) 
We note that |5]) is the simplest form of (\> a and <f> e providing dimensional homogeneity. 

Equation Q is often used in modeling viscoelastic bodies. It is used in [TT] for the study of 
the wave propagation in an unbounded domain. Equation ([7]) is used in [2, 13J as well as in [3[7], 
where the wave motion, stress relaxation and creep, are studied on a bounded domain. We also 
mention that the wave motion in a body, described by a more general model than ([5]) is studied 
in [T2] • We refer to O H3 for the detailed account of applications of fractional calculus in 
viscoelasticity. Problems similar to (JXJ) - (HJ were also treated in [T51 HH] with the constitutive 
equations related to the distributed-order model ([2]) in the special cases. 

We treated in [7] the creep test of a material described by the constitutive equation ([7]) 
and concluded (according to numerical examples) that the material is solid-like, while in [5] the 
constitutive function was given by 

(l + a D?)a{x,t) =E(b 00 rf° + 61 O of 1 + b 20 B^ e (x,t) , (9) 

where a, 60, b±, 62 are positive constants, < a < j3 < fi-y < j3 2 < 1, and the conclusion was 
(again based upon numerical examples) that the material is fluid-like. The constitutive equation 
([9"|) is proposed in |22) . In this work we treat numerically the creep test for solid- like materials 
described by 

Remark 1 Constitutive functions (or distributions) CT and <f> e in 0) have to satisfy the re- 
strictions following from the Second Law of Thermodynamics, see J7J [5]/. We refer to for 
a systematic review of restrictions if <p a and 4> e are given in the form of sums of the Dirac 5 
distributions 

N M 

$<J W := an 5 ~ ' ^Al) -=^Z b m5(l ~ P m ) 1 , Pm € [0, 1] . 

n— m— 
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Remark 2 Differences between solid and fluid-like materials are observed in the creep test (i.e. 
when a material is subjected to a sudden, but later constant force on its free end). Namely, solid- 
like materials creep to a finite displacement, while the fluid-like materials creep to an infinite 
displacement. 

The paper is organized as follows. In §[2] we write the system ([1]) - (J4J) in the dimensionless 
form and obtain (jTUJ) - (|13l) . Then we formally apply the Laplace transform to (JTUJ) - (flU]) , define 
the function M and obtain the solutions to (fTTJj) - (fT5|) in the Laplace domain via the forcing term 
and solution kernel. Section [3] is devoted to the verification that the function M in the cases 
of the fractional Zcner (0 and distributed-order model ((?]) satisfies assumptions, cited from [8], 
that imply the existence and uniqueness of the solutions to (fT0|) - (|13p . Then, we write theorems 
on existence and uniqueness of solutions, that are proven in [8]. The explicit form of the solution 
u, given in Theorem 21 is used in §[5] in order to plot the solution. The plots are given and 
discussed for the fractional Zener model and for two different forcing functions. 

2 Formal solutions 

We start from the system JTJ) - ([4]) and write it in the dimensionless form. Then, by the Laplace 
transform method, we obtain the displacement u and the stress a as the convolution of the 
external force F and solution kernels P and Q, respectively. Determination of P and Q will be 
given in §[3J 

The system QJ - Q transforms into 

d d 2 d 

—a(x,t) = n 2 —u(x,t), e(x,t) = — u(x,t) , x € [0, 1] , t > 0, (10) 

i ,i 

<M7)oDMM)d 7 = / <t> e (l)o^le(x,t)dj, x€ [0,1], * > 0, (11) 
Jo 

u(x,0) = 0, —u(x,0) = 0, a(x,0) = 0, e(x,0) = 0, xe[0,l], (12) 

d 2 

u(0,t) = 0, -<r(l,t)+F(t) = — «(l,t), t>0. (13) 
This is done by introducing the square root of the ratio between the masses of a rod and a body 



pAL 



and dimensionless quantities 



x ~ * - u _ a - (j>„ _ - _ 4> e p _ F 

.T<Pe- , rr\7- * ~ AE - 




In writing (1101) - (1131) we omitted bar over dimensionless quantities. Note that the choice of 
dimensionless quantities implies that the case of a rod without the attached mass (m = 0) 
cannot be studied as a special case of equations (fTTJ|) - (TTSj) . 

In order to solve the system (TlTJ|) . ((TTJ) subjected to the initial (j!2l) and boundary data (TT3l . 
we use the Laplace transform method. Recall that the Laplace transform of / € L\ oc (R) , / = 
in (—oo,0] and \ f (t)\ < ce kt , t > 0, for some k > 0, is defined by 

/>oo 

/ (s) =C[f (t)} (s) := / / (t) e- st dt, Res>k 
Jo 
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and analytically continued into the appropriate domain D. Moreover, we consider our problems 
within the the space of tempered generalized functions supported by [0, oo) , denoted by S',. The 
Laplace transform within this space is an extension of the classical one, given above. Namely, 
any g e S' + of the form g :— oD]f, 7 € (0,1), where / is as above (polynomially bounded) 
satisfies C [g {€)} (s) = s 7 / (s) , Re s > 0. We refer to [33] for the properties of elements of S' + and 
their Laplace transforms. 

Applying formally the Laplace transform to (|TU1) - (IT51) we obtain 

d d 

—a (x, s) — k 2 s 2 u (x, s) , e (x, s) = —u (x, s) , x £ [0, 1] , s £ D, (14) 
ox ox 



&(x,s) 4*a (7) « 7 d7 = e {x, s) / <j) e (7) s 7 d7, x £ [0, 1] , s £ D, 
Jo Jo 



(15) 



2(0,s) = 0, o- (1, s) + s 2 M (1, s) =F(s), sell. (16) 

By ([T5"[) we have 



gfos) = M^) g ^' s )' seD, (17) 



where we introduced 



M (s) := 



fn ^ (7) s7d 7 

, seD. (18) 



\ Jo <^ (7) S 7 d 7 

Thus, using ([5]) and ©, in the cases of constitutive equations ([5]) and ([7]) we obtain 



• U( - s) = «€C\(-oo,0], 0<a<6, a€(0,l), (19) 



M{8) = / ln ( bs ) as 1 S £C\ (-00, 0] , < a < b. (20) 
y In (as) bs — 1 

So, in the sequel D = C\ (— 00, 0] . 

In order to obtain the displacement u, we use (Tl4|) . (fl7|) and obtain 

<9 2 

— 2 (x, s) - k 2 s 2 M 2 (s) 2 (x, s) = 0, x G [0,1], s G C\(-oo,0], (21) 

The solution of ([2TJ is 

2 (x, s) = Ci (s) e KXsM M + C 2 (s) e - Ka:sM ( s ), x G [0, 1] , s G C\ (-00, 0] ; 
C\ and C2 are arbitrary functions which are determined from (fT6")h as 2C = Ci = — C2. Therefore, 
2 (x, s) = C (s) sinh (kxsM (s)) , x G [0, 1] , s £ C\ (—00, 0] . (22) 
By COO) 2, (HTJ) and (J22) we have 

cr(x,s) = C(s) j^fhj cosh (kxsM(s)), x G [0, 1] , s G C\ (-00, 0] . (23) 



Using (|22p and ([23)1 at x = 1, as well as (1161) 9 we obtain 

M (s) F (s) 

° ^ = s (sM (s) sinh (ksM (s)) + k cosh (ksM (a))) ' seC \(~°°'°]' 
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Therefore, the Laplace transforms of the displacement and stress from (|22l) and (|2"3"|) are 

u 0, s) = F (s) P (ar, s) , a (x,s) = F (s) Q (x,s) , x € [0, 1] , s e C\ (-oo, 0] . (24) 

where 

p. / 1 M (s) sinh (kxsM (s)) r „ „. . „. ,„ N 

P x ' s = ~ DTT^ITT » f r It mi W x£ °' X ' «€C\(-oo,0 , 25 

s sAJ (s) sinh (ksA? (s)) + k cosh (ksA4 (sJJ 

~ . . k cosh ( kxsM (s)) .„ „., . „. - „. 

Q(z,s) = % ■ w — ^^^7 — ttttt> xe 0,1, seC\ -oo,0 . (26 

v ; sAf (s) sinh (ksM (s)) + k cosh (ksM (s)) l j vv j v ; 

Applying the inverse Laplace transform to (|24|) we obtain the displacement and stress as 

it (x, t) — F (t) * P (x, t) , cr (x,t) = * Q (x,t) , z e [0,1], t> 0. (27) 

The validity of these formal expressions will be proved in the sequel. 

3 Explicit form of solutions 

In order to obtain the displacement u and stress a by (|27p . we have to obtain functions P and 
Q, i.e., to invert the Laplace transform in (|25p and (|2"r?)) . First, we examine the behavior of the 
function M, given by (fl8|) . in the limiting cases as |s| — > oo and \s\ — > (by this we mean only 
those s that belong to C\ (— oo, 0]) in the special cases when M takes any of the forms given by 
flUD and (HOD- 

If M is given by (HU) or ([20) we have 

|M(s)|«y|, as |s|->oo, and |M(s)|«l, as |s| -> 0. (28) 

We use results from [8] in order to obtain the displacement u and the stress a. In order to do 
so, we recall assumptions on M that have to be satisfied. 
We shall analyze a function of complex variable 

/ (s) := sM (s) sinh (ksM (s)) + n cosh (ksM (s)) , s € C. (29) 

Let M be of the form M (s) = r (s) + ih (s) , as \s\ — > oo. Then 

(Al) 

lim r (s) = Cqo > 0, lim h (s) = 0, lim M (s) = Co, 

|s| — >oo |s|— >oo s|— fO 

for some constants Cqo, cq > 0. 

Let s„ = + if n , n € N, satisfy the equation 

f(s)=0, s&V, (30) 

where / is given by (|2"9")l . Then: 

(A2) There exists no > 0, such that for n > no 

Ims„ eM + =>/i (s„) < 0, Ims„ e M_ =>■ /i (s n ) > 0, 
where /i := ImAf. 
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(A3) There exist sq > and c > such that 



-(sM {s) ) 



> c, \s\ > s . 



(A4) For every 7 > there exists 8 > and so such that 

|(s + As) M(s + As) - sM{s)\ < 7, if |As| < 6 and |s| > s . 

In order to write M = r + ih, as required above, we start from M 2 = u + iv and obtain the 
system 

r 2 — h 2 = u, 2rh = v. 
Solutions of the previous system belonging to the set of real numbers are 



V2 



r = ± — y \J u 2 + v 2 + u, 



± 



V2 



\J \/ u 2 + v 2 



Assume — — >■ 0, u > 0. Then by using the approximation (1 + a;) 
([3Tj) and ([32]), we have 

1 M 



(31) 
(32) 

1 + ax, as x — > + , from 



2a/u 



(33) 



Proposition 3 Functions M given by U9p and i20\) satisfy (Al) - (A4) . 
Proof. Consider M given by (|19p . If we write s = re ll<3 , by p^|) . we have 
1 + as Q 



M 2 (s) 



, a < b 
1 + 6s Q 

1 + (a + 6) r° cos (ay?) + afrr 2a 

1 + 26r a cos (a<^) + (br a ) 2 
a .b — a sin (cap) 



(b — a) r a sin (atp) 
1 + 26r a cos (aip) + {br a f 



b l ~ b 2 



as s — > 00. 



Using ([3"3"j) we obtain 



V 2y a b z r a 

Let ^ e (0, 7r), then sin (o^) > 0. This implies Re (M 2 (s)) > and Im (M 2 (s)) < 0. Therefore, 
we also have Re(M(s)) > and Im (M(s)) < 0. Similar arguments are valid if ip £ (— 7r,0). 
Hence, we finally have 



r(s) 



, h(s) 



1 b b — a sin («</?) 



as s — >• 00. 



6' 2 V a b 2 r a 

Next we prove that M in ([191) satisfies (Al) . By ([34]) and ([19]), we have 



(34) 



lim r (s) 

I s| — ¥00 



, . lim h(s) = 0, lim \M (s)\ = 1. 

|s|— >oo I s I — J-0 
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Validity of assumption (A2) follows from 

In order to show that M in (fT9)) satisfies (A3) , we use (fT9)) and obtain 



ds 



Thus, by (HP 



(sM (s)) = M (s) 1 - 



- (S M( S) ) 



a (6 — a) s Q 



2(I + as Q )(l + 6s a ) 



— , as |s| — > oo. 



(35) 



We have that there exists £ such that 

\(s + As) M (s + As) -sM{s)\ < \ As\ 



-(sM(s)) 



Since ^ (sM (s)) , by (1531) . i s bounded as |s| — > oo and if |As| < 6 for some > 0, we have that 
(A4) is satisfied. 

Now consider M given by (j20)) . With s = re' v we have 



, , 9 , . In (bs) as — 1 
M (s) = , ; '- -, a < b 



In (as) bs — 1 

In (ar) In (br) + ip 2 abr 2 — (a + b) r cos tp + 1 

In 2 (ar) + <(0 2 (&r) 2 — 26r cos (/? + 1 In 2 (ar) + </? 2 (for) 2 — 26r cos ip + 1 



(6 — a) r siny; 



. / In (ar) In (for) + </? 2 (b — a) r sin 



</5 In | aor 2 - (a + 6) r cos </? + 1 



In (ar) + ip 2 (br) —2brcos<p + l In (ar) + p 2 (br) —2brcosip + l 



a a b ip 
- — l— In ^ — 



& a In (ar) ' 
Using (j3"3"j) we obtain 



as s 



oo. 



"2 V 6 a In 2 (ar) 

Let S (0, tt), then sin (a<p) > 0. This implies Re (M 2 (s)) > and Im (M 2 (s)) < 0. Therefore, 
we also have Re(M(s)) > and Im (M (s)) < 0. Similar arguments are valid if ip € (— 7r,0). 
Hence, we finally have 



/a 1 [a b p 

r(s) = W— , n(s) = — — \l — In 



6 ' ft ^ 2 V 6 *" a In 2 (ar) ' 



(36) 



Next we prove that (J20J) satisfies (Al) . By (J36J) and ([20]), we have 



lim r (s) 



lim h(s) = 0, lim |M (s)\ = 1. 

b \s\— >oo I s I — J-0 



Validity of assumption (A2) follows from ((3"6"]) . 

In order to show M in ((2T)]) satisfies (A3) , we use (j2"0")) and obtain 



-(sM( S )) = M( S )(l- 



hv 



(b — a) s 



2 In (as) In (6s) 2 (as — 1) (bs — 1) 
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Thus, by (PI) 



(sM(s)) 



as s — > oo. 



(37) 



Using the same arguments as above, we have that (A4) is satisfied, since sM (s) , s € V, by 
(|3T)l . has bounded first derivative. ■ 

The existence and the uniqueness of u and cr, as solutions to system (|10[) - (|13p is guaranteed 
by the fact that M in all two cases satisfies (Al) - (A4) . Recall, / is given by (j!29"]) and s n , n e N, 
are solutions of (I3U1) . The multiplicity of zeros s n is one for n large enough. 

Theorem 4 ( [8] ) Let F £ S' + and suppose that M satisfies assumptions (Al) - (A4) . Then the 
unique solution u to H10\) - \13\) is given by 



where 

P{x,t) 



Im 



u (x, t)=F{t)*P (x, t), x e [0, 1] , t > 0, 
M (qe- i7T ) sinh (nxqM (qc^)) 



-qi 



qM (qe~ i7r ) sinh {nqM (qe~™)) + k cosh (nqM (ge~ i7r )) J q 

OO 

2 Re ( Res I? ( x ' s ) eS *' s «) ) > xe[o,i], t> o, 



-dq 



P(x,t) 
The residues are given by 



n=l 

0, x£ [0,1], t<0. 



Res 



(P (x, s )e st ,s„) 



1 M (s) sinh (kxsM (s)) 



x e [o, l] , t > o, 



T/ien P £ C ([0, 1] x [0, oo)) and u £ C ([0, 1] , 5^) . In particular, if F £ L\ oc ([0, oo)) , then 
u is continuous on [0, 1] x [0, oo) . 

The following theorem is related to stress a. We formulate this theorem with F = H, where 
H denotes the Heaviside function, while the more general cases of F are discussed in Remark [51 
below. 

Theorem 5 (|8j) Let F — H and suppose that M satisfies assumptions (Al) - (A4) . Then the 
unique solution an to M0\) - H3\) , is given by 



a H (x,t) = H(t) + -f Im(- 



cosh i^nxqM (qe 17T )) 



-qi 



qM (qe in ) sinh (nqM (ge i7T )) + Kcosh {nqM (qe i7T )) 

+2^Rc(Rcs(ct h (x,s)e st ,s n )) , x £ [0, 1] , t > 0, 
n=l 

a H (x,t) = 0, x£ [0,1], t< 0. 
The residues are given by 



-dq 



(38) 



Res (o~ u (x,s)e st ,s n ) 



k cosh (kxsM (s)) st 



, ie[o,i],{>o. 



In particular, o~h is continuous on [0, 1] x [0, oo) 
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Remark 6 Q8J) 



1. The assumption F = H in Theorem^ can be relaxed by requiring that F is locally integrable 
and 

for some a £ (0, 1) . This condition ensures the convergence of the series in i38\ ). 

2. If F = 5, or even F (t) = 4r%8 (t) , one uses an, given by i38\) . in order to obtain a as the 
k + 1-th distributional derivative: 

d k+i 

4 Numerical examples 

The displacement u as a solution to (|10p - (IT51) is given in Theorem 2] We present various 
numerical examples for constitutive models: fractional Zener and distributed-order model of a 
solid-like viscoelastic body that are distinguished by the form of M: ([T5]l and PD]l. respectively. 

4.1 The case F = 5 

In order to plot time dependence of the displacement u for the several points of the rod as well 
as for the body attached to its free end, we chose the fractional Zener model and the force acting 
on the body to be the Dirac delta distribution, i.e. F = 6. We fix the parameters describing 
the rod: a = 0.2, b = 0.6, a = 0.45 and also fix the ratio between the masses of rod and body 
h := k 2 = 1. Plot of u as a function of time t for various points of a rod is shown in Figure [5J It is 

u(x,f) 




Figure 2: Displacement u(x,t) in the case F = S for x = 1 as a function of time t at x € 
{0.25,0.5,0.75,1} for t G (0,40). 

evident that the oscillations of the rod and a body are damped, since the material is viscoelastic. 



10 



One notices that initially there is a transitional regime of the oscillations. Afterwards, the curves 
resemble the curves of the damped linear oscillator. 

In order to examine the transitional regime more closely, in Figures [3] -E] we present the plots 
of u for smaller values of time, but for different values of x £ {0.5, 1,2}. We notice that the 

«(*>') u(x,t) 




Figure 3: Displacement u(x,t) in the case F = Figure 4: Displacement u(x 7 t) in the case F = 
5 for x = 0.5 as a function of time t at x £ S for x = 1 as a function of time t at x £ 
{0.25, 0.5, 0.75, 1} for t £ (0, 9). {0.25, 0.5, 0.75, 1} for * £ (0, 10). 

u{x,t) 




Figure 5: Displacement u{x, t) in the case F — 
5 for x = 2 as a function of time t at x £ 
{0.25,0.5,0.75,1} for t £ (0,11). 



shape of the curves depends on the ratio between the masses x, while later, the shape resembles 
to the shape of curves for damped oscillations. It could be noticed that regardless of the value 
of x there is a delay in starting oscillation for the points that are further away from the free end 
of a rod. This is due to the finite speed of wave propagation. Namely, the body [x = 1), which 
is subject to the action of the force, starts oscillating at t = 0, while the delay in the starting 
time-instant of the oscillation is greater as the point is further from the point where the force 
acts. Moreover, we see that different points of the rod do not come to their initial position, for 
the first time, at the same time- instant (see Figure [5]). This, as well as the initial delay depend 
on the mass ratio x. For the influence of x on the initial delay compare Figures - EJ Later on, 
again depending on x, the motion of the points become synchronized. 

Figures |6] and [7] present the plots of displacement u for fixed point of the rod x = 0.5 if 
the ratio between masses x varies. One sees from Figure [5] that the larger the mass of a rod is 
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Figure 6: Displacement u(x, t) in the case F — 8 at x — 0.5 as a function of time t € (0, 40) for 
x = 0.5 - dot-dashed line, x = 1 - dashed line and x = 2 - solid line. 

(then the value of >r is greater) the quasi-period (time between two consecutive passage of a fixed 
point through its initial position) of the oscillations is greater. That is due to the increased rod's 
inertia. Also, for larger times there is no significant influence of x on the heights and widths of 
the peaks which indicates that the damping effects are only due to the parameters a, b and a 
figuring in the constitutive equation. Figure [7] shows that the shape of the curve in transitional 
regime strongly depends on x. Moreover, the delay in the oscillations increases as the x increases 
which indicates that the speed of the wave propagation depends on the mass ratio. 



The aim of this section is the qualitative analysis of the behavior of a displacement u when there 
is a force, given in the form of the Heavisidc function, i.e., F = H, acting at the attached body. 
Thus, our results correspond to a creep experiment. The rod is modelled by the fractional Zener 
model, i.e., the function M is given by (|19l) . The parameters describing the rod are: a = 0.1, 
b = 0.9, a = 0.8. We present plots of u for the ratio between the masses of rod and body x=l. 

Figure [5] present the long-time behavior of displacement u. One notices that the rod creeps 
to a finite value of displacement so that lim^oo u (x, t) = x, x £ [0,1]. Figure [5] present the 
short-time behavior of u. We see that there is a delay in starting time-instant of a point of a 
rod. 

Figures [TU1 and [TT1 present the plots of time evolution of displacement u of a rod described by 
the Zener model for fixed point of the rod x = 0.5 if the ratio between masses x varies. Here, the 
parameters are a = 0.2, b = 0.6, a = 0.45. One sees, Figure [TO] that there is no dependence of 
the finite value of displacement in creep on the value of the mass ratio x. We see, Figure [TTJ that 
for small times there is an influence of x on the height of the peaks such that its height increases 
as x increases. This is due to the inertia, while for larger times the viscoelastic properties of the 
rod prevail. Similarly as in the previous section the delay in the oscillations starting time-instant 
increases as x increases. 
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4.2 The case F = H 
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Figure 7: Displacement u(x, t) in the case F — 8 at x — 0.5 as a function of time t € (0, 11) for 
x = 0.5 - dot-dashed line, x = 1 - dashed line and k = 2 - solid line. 
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Figure 10: Displacement u(x,t) in the creep experiment at x = 0.5 as a function of time t € 
(0, 100) for h — 0.5 - dot-dashed line, x = 1 - dashed line and k = 2- solid line. 
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